Testing the Warm Dark Matter paradigm with large-scale structures 
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We explore the impact of a AWDM cosmological scenario on the clustering properties large-scale 
structure in the Universe. We do this by extending the halo model. The new development is that 
we consider two components to the mass density: one arising from mass in collapsed haloes, and the 
second from a smooth component of uncollapsed mass. Assuming that the nonlinear clustering of 
dark matter haloes can be understood, then from conservation arguments one can precisely calculate 
the clustering properties of the smooth component and its cross-correlation with haloes. We then 
explore how the three main ingredients of the halo calculations, the halo mass function, bias and 
density profiles are affected by WDM. We show that, relative to CDM, the halo mass function is 
suppressed by 50%, for masses ~ 100 times the free-streaming mass-scale Mf B . Consequently, the 
bias of low mass haloes can be boosted by as much as ~ 20% for 0.25 keV WDM particles. Core 
densities of haloes will also be suppressed relative to the CDM case. We also examine the impact 
of relic thermal velocities on the density profiles, and find that these effects are constrained to 
scales r < 1 /i _1 kpc, and hence of little importance for dark matter tests, owing to uncertainties 
in the baryonic physics. We use our modified halo model to calculate the non-linear matter power 
spectrum, and find that there is significant small-scale power in the model. However, relative to the 
CDM case the power is suppressed. The amount of suppression depends on the mass of the WDM 
particle, but can be of order 10% at k ~ 1 h Mpc -1 for particles of mass 0.25 keV. We then calculate 
the expected signal and noise that our set of AWDM models would give for a future weak lensing 
mission. We show that the models should in principle be separable at high significance. Finally, 
using the Fisher matrix formalism we forecast the limit on the WDM particle mass for a future 
full-sky weak lensing mission like Euclid or LSST. With Planck priors and using only multipoles 
I < 5000, we find that a lower limit of 2.6 keV should be easily achievable. 
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I. INTRODUCTION 

Within the last two decades cosmology has experienced 
a major observational revolution. The measurement of 
the Cosmic Microwave Background (CMB) anisotropics, 
through such surveys as COBE and WMAP, and the 
mapping of the large-scale structures (LSS), through 
galaxy redshift surveys like 2dFGRS and SDSS, have had 
a great impact on our knowledge of the Universe. In par- 
ticular, it has led to the establishment of the 'standard 
model' for cosmology: ACDM model [l[ - meaning a 
Universe whose mass density is dominated by collision- 
less Cold Dark Matter (CDM) relic particles, but whose 
total energy density today is dominated by the constant 
energy of the vacuum A, which is driving the current 
accelerated expansion of the Universe. However, funda- 
mental questions remain unanswered: What is the true 
physical nature of the cold dark matter? and what is the 
true origin of the accelerated expansion? In this paper we 
focus on developing tests for the dark matter model based 
upon the influence that different particle candidates have 
upon the statistical properties of the LSSs. 

There is a large body of indirect astrophysical evidence 
that strongly supports the CDM hypothesis - i.e. a 
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heavy, cold thermal relic particle, e.g. weakly interacting 
massive particles, etc., that decoupled from normal mat- 
ter very early in the history of the Universe. However, de- 
spite the many successes of ACDM, there is also a grow- 
ing body of evidence in conflict with parts of the hypoth- 
esis. Firstly, hi-resolution simulations of CDM predict 
several hundred dwarf galaxy halo satellites within the 
Local Group, whereas the observed abundance of dwarf 
galaxies is a factor of 10 lower although counting 

and weighing local group galaxies is not easy [l(| [nj. 
Secondly, the inner density profiles of haloes in CDM 
simulations are cuspy, with the inner density falling as 
p(r) oc r _1 [1343 1 whereas the density profiles inferred 
from galaxy rotation curves are significa ntly shallower 
[HI (and for more recent studies see [l7l42C| and refer- 
ences there in). Thirdly the observed number of dwarf 
galaxies in the voids is far smaller than expected from 
CDM 0. 

Over the past decade much work has been invested 
in attempting to understand these discrepancies. Some 
of these problems may be reconciled through more care- 
ful treatment of baryonic physics and galaxy formation. 
More exotic explanations for some of these problems 
would require modifications to gravity or scalar fields 
that operate only in the dark matter sector 0, HH . Per- 
haps a less dramatic explanation might be the modifica- 
tion of the physical properties of the dark matter particle. 

Over the past decade, several modifications to the 



CDM hypothesis have been forwarded [see for exam- 
ple [22| . In this work we focus on Warm Dark Matter 
(WDM) in a low-density A dominated universe (hereafter 
AWDM). In this scenario, the dark particle is considered 
to be lighter than its CDM counterpart, and retains some 
stochastic primordial thermal velocity, whose momentum 
distribution obeys Fermi-Dirac statistics |23l. [23 |. Two 
potential candidates are the gravitino [25j and a ster- 
ile neutrino (2(| , both of which require extensions of the 
standard model of particle physics [23|, [27j[. The main 
cosmological effect of WDM is that, at early times while 
it remains relativistic, density fluctuations are suppressed 
on scales of order the Hubble horizon at that time - and 
particles are said to free-stream out of perturbations. 
In fact, until matter-radiation equality, the WDM par- 
ticles can also still travel a comparable distance whilst 
non-relativistic. Consequently, the matter power spec- 
trum becomes suppressed on small scales. At later times 
the particles cool with the expansion and at the present 
time essentially behave like CDM - modulo any remain- 
ing primordial velocity dispersion [for more details see 

Past research on the impact of the WDM hypothesis 
on the phenomenology of cosmic structures has primarily 
focused on small-scale phenomena [id |28- 32] , although 
there are some notable studies of large-scale structures 
[23l [33|. In short, these results show that halo abun- 
dances for masses below the free-streaming mass-scale 
are suppressed relative to CDM. The amount of suppres- 
sion is still a point of debate: initial studies, suggested 
that low mass haloes might form through the fragmenta- 
tion of larger perturbations, in particular along filaments 
23]. However, recent numerical studies f34ll35j| have cast 
doubt on this as a viable mechanism, since the abundance 
of structure forming below the free-streaming mass-scale 
appears to be correlated with the initial particle grid and 
inter-particle separation. Some simulation results have 
also claimed that the density profiles of WDM haloes are 
shallower in the inner core regions [Hj]. Further, whilst 
the linear matter power spectrum is suppressed, the non- 
linear one regenerates a high fc-tail |29j |. 

Observational constraints suggest that sterile neutri- 
nos can not be the dark matter: the Lyman alpha for- 
est bounds are m Vs > 8keV [13, [36[ , whilst those from 
the X-ray background are m Ue < 4keV 37|. However a 
more recent assessment has suggested that a better mo- 
tivated particle physics model for the production of the 
sterile neutrino may evade these constraints, with the 
new bound m Ve > 2keV [Hj]. It therefore seems that 
independent methods for constraining the properties of 
WDM from cosmology would be valuable. 

In a recent paper, [33, it was proposed that the 



paper, 

AWDM scenario could be tested through weak lensing 
by LSS. The advantage of such a probe is that it is only 
sensitive to the total mass distribution along the line of 
sight. However, to obtain constraints on the WDM par- 
ticle mass, an accurate model for the nonlinear matter 
clustering is required. Markovic et al. [39[ attempted to 
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FIG. 1: Free-streaming mass-scale (Mf s ) as a function of the 
mass of the WDM particle mass towdm. Haloes that have 
masses M > Mf B may form hierarchically (hatched region). 
Haloes with M < M{ a (empty region) will have their peaks in 
the primordial density field erased and so may not form hier- 
archically. Haloes in this region may possibly form through 
fragmentation. The solid red band shows the mwDM allowed 
by the Lyman alpha forest [2tJ (note that we have rescaled 
m V3 — > mwDM using Eq. ©). 



construct this in two ways, firstly through use of the halo 
model and secondly through adapting a fitting formula 
for CDM power spectra. There appeared to be some dis- 
crepancy between the predictions of the two approaches. 
Furthermore, their halo model calculation suggested that 
the bias of high mass haloes was significantly larger than 
for the case of CDM. One of the aims of this paper is to 
construct an improved halo model calculation for appli- 
cation to this problem. 

The paper is broken down as follows. In <JTT]we briefly 
review the impact of WDM on the linear theory evolu- 
tion of structure. In ^Illl we describe our new halo model 
approach. In S jlVI we discuss the ingredients of the halo 
model in the context of WDM. In fJV]we show our results 
for the clustering statistics and the weak lensing observ- 
ables. In j}V CI we forecast a limit on the particle mass, 
obtainable from a future full-sky weak lensing survey. In 
Q VII we summarize our findings and conclude. 



II. THEORETICAL BACKGROUND 

One of the major differences between structure forma- 
tion in the ACDM and AWDM scenarios, is that whilst 
the dark matter particle is relativistic it free-streams 
through the Universe and so suppresses the growth of 
density perturbations on scales smaller than the free- 
streaming length [H, H, S3, EJ [73]. The lighter the 
WDM particle the longer it remains relativistic and hence 
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the larger the free-streaming scale. Following conven- 
tion, we shall define the free-streaming scale as the half- 
wavelength of the mode for which the linear per turbation 
amplitude is suppressed by a factor of 2 [23l. l30t |4~D ] . Fit- 
ting to the solutions of integrations of the collisionless 
Einstcin-Boltzmann equations through the epoch of re- 
combination, sets the free-streaming scale to be [42j]: 
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where fiwDM is the total energy density in WDM, rel- 
ative to the critical density for collapse, and TOwdm is 
the mass of the dark matter particle in keV. From the 
scale Af s we may infer a putative halo mass whose for- 
mation is suppressed by the free-streaming. In the initial 
conditions this is given by [30l |: 
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In Fig. [I] we show the relation between the free- 
streaming mass-scale and the mass of the WDM particle 
candidate for our adopted cosmological model. Haloes 
that would have masses below the free-streaming mass- 
scale will have their peaks in the primordial density field 
erased. The impact of the free-streaming of the WDM 
particles on the linear growth of structure can be further 
described by the matter transfer function. This can be 
represented (23j, [Hj], 



FIG. 2: Linear mass power spectra as a function of 
wavenumber in the WDM and CDM scenarios. The solid 
black line shows results for CDM. The dot-dashed lines 
represent the WDM power spectra for the cases where 
rawoM = m x e {0.25, 0.5, 0.75, 1.0, 1.25} keV. The lighter 
the WDM particle the more the power is damped on small 
scales. 



III. HALO MODEL APPROACH 



A. Statistical description of the density field 
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where -Pl^ dm and Pfj^ M are the linear mass power spec- 
tra in the WDM and CDM models respectively. Note 
that in order to avoid cumbersome notation, we shall 
adopt the convention Pu n = ^Lhi DM ' likewise Pnl = 
Pnl >M f° r the nonlinear power spectrum in the WDM 
model. Fitting to results of full Einstein-Boltzmann in- 
tegrations gives p = 1.12, and a of the order [4lj |: 
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In Fig. [5] we show how variations in the mass of the 
WDM particles effects the linear matter power spectrum. 
Clearly, the lighter the WDM particle the more the power 
is damped on small scales. 

Note that in the above we are assuming that the WDM 
particle is fully thermalized, i.e. the gravitino. Follow- 
ing [4l|, this can be related to the mass for the sterile 
neutrino through the fitting formula: 



4.43keV 



ra WDM \ 4/3 / w wdm \ ~ x / 3 
IkeV / V 0.1225/ 



(5) 



We would like to know how the nonlinear matter power 
spectrum changes for WDM models. In order to do this 
we shall adapt the halo model approach. For the case of 
CDM, the halo model would assume that all of the mass 
in the Universe is in the form of dark matter haloes and 
so the density field can be expressed as a sum over haloes. 
In WDM models this ansatz can be no longer considered 
true, owing to the suppression of halo formation on mass- 
scales smaller than Mf s . We modify this in the following 
way: we suppose that there is some fraction of mass in 
collapsed dark matter haloes and that there is a corre- 
sponding fraction in some smooth component. Hence we 
write the total density of matter at a given point as: 



p(x) = p s (x) 
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where p s is the density in smooth matter, and U(x, Mi) = 
p{x\M)/M is the mass normalized isotropic density pro- 
file, and both depend on the mass of the WDM particle 
- for convenience we suppress this label. The sum in 
the above extends over all N haloes and subhaloes. The 
mean mass density of the Universe can be written, 



(p(x)) =p = p s +ph 



(7) 
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where the expected total mass density in WDM haloes 
and the smooth component can be written, 



Ph = fp 



Ps = (1 - f)f> 



(8) 



where these terms have the explicit forms [431 ] : 
&r) = - / dMM 2 n{M) 

Ph JM mt 



where p is the mean density of all matter, and / the 
fraction of mass in haloes. This can be written, 



1 f°° 

f = - dMMn(M) 

P A/ cut 



(9) 



where n(M) is the number density of WDM haloes per 
unit mass and M cut is the mass-scale below which there 
are no haloes. In this study we will take M cut = Mf s , but 
will formally make the distinction throughout. 

Let us now work out the clustering statistics in this 
model. The two point correlation function of the mass 
density can be written: 

(p(x)p(x + r)) = ([ i o s (x)+p h (x)] [p s (x + r)+p h (x + r)]) 
= (p s (x)/9 s (x + r)) + (p s (x)p h (x + r)) 
+ (Ph(x)p s (x + r)) + (/> h (x)ph(x + r))(10) 

Owing to the statistical homogeneity and isotropy, we 
have that (p s (x)ph(x + r)) = (/9h(x)p s (x + r)). The cor- 
relation function between two random density fields X 
and Y can be written: 

6c Y (|r|) = <(5x(x)Mx + r)) , (11) 
where the dimensionless over-density field S is defined: 

Px - Px 



Sx = 



Px 



(12) 



On using the above definition, we find that 
fair) = (l-ff^(r) + 2f(l-f)U(r) + f 2 ^ h (r) . (13) 

Thus the correlation function is simply the weighted av- 
erage of the smooth-mass correlation function £b, the 
mass weighted halo-halo correlation functions £hh and the 
smooth mass-halo cross correlation function £hs- In or- 
der to proceed further, we need to understand these three 
correlation functions. 

B. Halo- halo correlation function 

Let us first consider the halo-halo correlation function, 
this can be derived following [43l and references therein] . 
The clustering is the sum of two terms: the 1-Halo term, 
which describes the clustering of points within the same 
halo and the 2-Halo term, which correlates the mass dis- 
tribution between haloes. Thus we have, 



xj d 3 x£/(x|Af)[/(|r-x|,M) ; (15) 

pi t\ LA/ cut 
x J d a mU(Ti - m\M)\ £; e h nt (|xi - x 2 |) , (16) 

where in Eq. dTp]) xi and x 2 denote the position vectors 
of the center of mass of halo 1 and 2, and ri and r 2 
denote position vectors of the points to be correlated, 
respectively. Thus ri — xi gives the radial vector from 
the center of halo 1 to the point 1, and likewise for halo 2. 
Further, owing to the statistical isotropy of the clustering 
£(r)=£_(|r|)=£(|ri-r 2 |)._ 

In principle the correlation function of halo centers 
^ lt {y\M 1 ,M 2 ) (with y = |xi — x 2 |), is a complicated 
scale dependent function of Mi, M 2 and y [e.g. see [44|. 
We assume that there is a local deterministic bias rela- 
tion between haloes and dark mass, smoothed on a scale 
R. Hence, we have lii-lil: 



h(M) 



[Mx)r , (17) 



where the bi(M) are the nonlinear bias parameters. Thus 
the halo center correlation function has the form: 

&(V) = 6i(Mi)6i(M 2 )e K (y) 

+1 [6i(Mi)6 3 (M 2 ) + 68(M 1 )6 1 (M 2 )] (M x )<^( x + v)) x 
^^^^)(4(x)4(x + y )) x + ...(18) 



^2H 



(14) 



where £r(v) is the nonlinear matter correlation function 
smoothed on scale R. On truncating the above expression 
at first order in bi and inserting into Eq. (|16p . one finds: 

1 2 f f°° 
&h(r) = — TT / dM t MMMi)n{Mi) 
PI t =1 Um c „ 

x J d^f/Cri-Xi^Ufldxi-xa])^) 

Before moving on, we note that in practice, in or- 
der to make accurate predictions in the halo model on 
small scales, we must take into account halo exclusion 
[lHIU 48]. That is, the separation of two distinct haloes 
can not be smaller than the sum of their virial radii. Oth- 
erwise, they would have been identified as a larger halo. 
However, when modelling WDM there is a new "exclu- 
sion" condition that arises, the smooth component is not 
really smooth. Formally there are holes of zero density 
in the smooth distribution where haloes are. Since these 
calculations are first order attempts, we shall leave such 
detailed analysis for future study. 
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C. Smooth-smooth correlation function 

Let us next consider the smooth mass component of 
the model. We shall assume that this represents all of 
the mass that was initially in mass peaks of the density 
field that were M < Mf s . In direct analogy with the den- 
sity field of dark matter haloes, we may assume that the 
smooth component can be related to the total mass den- 
sity field through a deterministic local mapping. Hence, 



Wx),j s [Mx)] = Ef[feW]' 



(20) 



where the b s> i are the nonlinear bias parameters of the 
smooth component. This leads us to write an expression 
similar to Eq. (TT8l) . 



&, R (r)=bl 1 Ur)+0(b 



s,2, U s ,3, 



D. Halo-smooth mass cross-correlation 



(21) 



The model also requires us to specify the cross- 
correlation between the mass in haloes and the mass in 
the smooth component. This can be written in the fol- 
lowing fashion, 

Chs(r) = — / dMMn(M) [ d 3 xit/(ri - Xi \M) 
Ph Jm cM J 

x6i(M)6 s , 1 ^(x 1 - r 2 ) + <D(b 2 , b s<2 ...) (22) 
where r = \r\ — r%\. 

E. Power spectrum 

The correlation functions are related to the matter 
power spectra through a Fourier transform [43|, [49|, [5(| : 



P(k) = J d 3 r£(r)e 



d 3 k 

(2tt) 3 



P(fc)e" 



(23) 
one finds 



On applying this to Eqns (TJJ}, (gT]), and 

Pss(k) = (l-/) 2 P fi (fc)+2(l-/)/P sh (A : )+/ 2 P hh (fc) (24) 

where the power spectra of the smooth-smooth, smooth- 
halo and halo-halo components, can be written to first 
order in the halo and smooth density bias parameters as, 



P ss — fr s ,i-PNL : R 
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(25) 

dMMb 1 {M)n{M)U{k\M){^) 



M , 
oo 

Af out 



dM Mbi(M)n(M)U (k\M) 



, , dMn(M)M I U 2 {k\M) 

Ph JMeut 



(27) 



where Pnl,r = -Pnl,r(^) is the nonlinear mass power 
spectrum smoothed on a scale R, and U(k\M) is the 
Fourier transformed mass normalized density profile, 
U(x,Mi). 



F. Bias of the smooth component 

Absent so far is any explanation as to how one might 
infer the bias factors of the smooth component b s ^. We 
shall now use a mass conservation argument to derive a 
relation between the bias of the smooth mass and the 
bias of the haloes that involves no new free parameters. 
That is, if we fully understand halo biasing, then we may 
immediately obtain a full understanding of the density 
distribution of the smooth density component. 

Consider again the density at a given point, this may 
be written as in Eq. ©, i.e. p(x) = p s (x) + jOh(x). On 
using Eqns ©, © and (TI"2")) . then we may write this in 
terms of the overdensities in the smooth and clumped 
matter as 



<5(x) = (l-/)5 s (x) + /5 h (x) 



(28) 



The mass density associated with matter in haloes can be 
written as an integral over all haloes above the mass-scale 
M cut in the following manner, 

/Oh(x) = Ph [1 + 4(x)] 

dMMn(M) [1 + <5 h (x|M)] . (29) 

M cut 

If we assume the local model for halo biasing of Eq. (|17[) , 
then the density field of haloes takes the form: 

&(x) = ^ !°° dMMn{M)Y^ b -^[8 R (*)] 1 (30) 

Ph JM cat V l - 

where ph = ff>- Next, let us define the effective mass 
weighted halo bias parameters as, 



-off 



J Mcut dMMn(M)bi{M) 



Then we may re-write Eq. (|30D more simply as 
4(x) = ^^Wx)f . 



(31) 



(32) 



Substituting the above equation and the local bias model 
from Eq. flU]) into Eq. (gHJ), we then find 



£[(l-/Xi + /&f] 



(33) 
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On inspecting the above equation, two conditions may 
be immediately noticed. Firstly, in order for the left and 
right sides of the equation to balance, we must have that 



1 = (1 - fib,,! + fb\ 



cff 



(34) 



This leads us to the result that the bias of the smooth 
component must be given by, 



&s,l = 



1 - fbf 

1-/ 



(35) 



Secondly, if the smooth component is indeed linearly bi- 
ased with respect to the total matter overdensity, then we 
must also have that the nonlinear bias coefficients vanish 
when integrated over all haloes: 



dMMn(M)bi(M) = (i ^ 1) 



(36) 



On the other hand, supposing that the smooth compo- 
nent is not simply linearly biased, but is nonlinear ly bi- 
ased a la the local model, then the above condition is 
relaxed to the condition: 



(f - f)b s ,i + fbf = (i ± 1) 



This leads us to the condition, 



b s .i = - 



i-/ 



(Mi) 



(37) 



(38) 



Hence, if the haloes and the smooth mass component are 
locally biased, then all of the nonlinear bias parameters 
for the smooth component can be derived from knowl- 
edge of the halo bias parameters. In ^IVBI we will show 
how the bias of the smooth component varies for different 
mass-scales of the WDM particle. 



where dn — n(M)dM is the abundance of WDM haloes 
of mass M in the interval dM, S c (z) = 1.686/Z3(z) is the 
linear theory collapse barrier for the spherical collapse 
model, and D(z) is the linear theory growth function. 
Note that we take the evolution of S c (z) and D(z) to be 
unchanged from their forms in the ACDM model. Our 
reasoning is that, if one follows the standard derivations 
for D(z) and S c (z), then, provided the WDM particles 
decouple early and are non-relativistic today, these func- 
tions are solely determined from the dynamics of the ex- 
pansion rate of the Universe. This assumption will of 
course break down at very early times, when the primor- 
dial thermal velocities are more significant and when the 
particles are relativistic and so growth suppressed. 

For the function f(v) we adopt the model of Sheth and 
Tormen 15211. and this has the form: 



2 



/stM = A(p)J-y/qv [1 + (yV)- 2p ] exp 



_ (40) 

For the ST mass function the amplitude parameter A is 
determined from the constraint § d\ogvf{v) = 1, which 
leads to A-^p) = [I + 2~*T [0.5 - p] /V [0.5]]. The pa- 
rameters are: {A = 0.3222, p = 0.3, q = 0.707}. The 
variance on mass-scale M and its derivative with respect 
to the mass scale can be computed from the following 
integrals: 



a 2 (M) 

d\oga 2 (R) 
d\ogM 



d 3 k 

(2tt) 3 
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Pun(k)W 
d 3 k 



3a 2 (i?) J (2tt) 



(41) 

Pun(k)WW . (42) 



where the filter functions are W = W(kR) and W = 
dW (kR) I 'dlog(kR) . For a 3D real-space spherical top- 
hat filter, transformed into Fourier space, we have: 



IV. HALO MODEL INGREDIENTS IN WDM 
MODELS 

In this section we detail the halo model ingredients and 
show how they change in the presence of our benchmark 
set of WDM particle masses. 



A. Halo mass function 

In order to predict the abundance of haloes in the 
WDM cosmological model, we shall adopt the Press- 
Schechter (PS) [51] approach, with some small modifi- 
cation. For mass-scales M > Mf s , we shall assume that 
the standard formalism holds. In this case we may calcu- 
late the abundance of objects through the usual recipe, 
only being sure to use the appropriate linear theory power 
spectrum for WDM. Hence, 



dn _\P_f( \ d lo S °" 2 
d log M ~ 2 M * ^ d log M 



S c (z) 
a(M) 



Wth (y) = o [sin y — y cos y] ; 

y A 



(43) 



(39) 



W^ H (y) - - [(y 2 - 3) sin y + 3y cosy] . (44) 

yA 

For mass-scales M < Mf s , free-streaming erases all 
peaks in the initial density field, and hence peak theory 
should tell us that there are no haloes below this mass 
scale. A blind application of the PS approach leads one 
to predict, suppressed, but significant numbers of haloes 
below the cut-off mass. 

In order to demonstrate this point explicitly, let us 
consider the following toy-model calculation. Suppose 
that the matter power spectrum is a simple power-law 
P(k) oc k n , and let us suppose that if the dark mat- 
ter particle is warm then the effect of the damping on 
the fluctuations due to free-streaming can be represented 
by a simple exponential damping. Hence, Avdm(^) = 
A s (k/k c ) n exp[— k/k c ], where k c is a characteristic damp- 
ing scale; perhaps k c = 27r/Af s . Let us now consider the 
behaviour of Eqs (gf|) and ((321) in the limit that R -> 0. 
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Mass [h _1 M 0] Halo Mass/M fs 

FIG. 3: Left panel: Halo mass function as a function of halo mass. Lines show the theoretical predictions from the ST model 
of Eq. (|40l) . The black solid line denotes the result for CDM; the dot-dash lines denote results for WDM for our standard set 
of particle masses (jtiwdm = rnx = {0.25, 0.5 0.751.0, 1.25} keV) with no cut-off mass applied; the solid colored lines denote 
the same, but with a cut-off mass-scale as given in Eq. ()49[) . and here we use cr\ og M = 0.5. Right panel: Fractional difference 
between the WDM and CDM mass functions, as a function of halo mass scaled in units of the WDM free-streaming mass-scale 
Mf s . Line styles are the same as in left panel. 
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FIG. 4: Left panel: Halo bias as a function of halo mass for various dark matter models. The solid black line denotes the 
results for CDM; the dot-dashed line denotes the halo bias for our standard set of WDM particle masses (jtiwdm = mx = 
{0.25, 0.5, 0.75, 1.0, 1.25} keV). The dashed lines denote the bias of the smooth component of mass as given by Eq. (|35[) . 
Right panel: fractional difference between the bias in WDM and CDM models as a function of halo mass-scaled in units of the 
free-streaming mass-scale Mf s . 



For the variance term we have: where we have used the fact that lim^^o W(kR) = 1. 

/d 3 k 
7^AvDM(fc) lim W 2 (kR) 

) 

^r[n + 3] (45) 



d 3 k 

AnA s kl 
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For the derivative of the variance we have: 



lim dl °f iR) 
R^o d log M 



lim — — — - 

R^o 3a 2 (R) 



d 3 k 



WDM 



lim W(kR) 

R-+0 



lim W'(kR) 

R-+0 



= -1 



(46) 

where we used the fact that liiriR^o W'(kR) = —3/2. 
Thus as M — > (or equivalently R 0), the mass func- 
tion of dark matter haloes becomes: 



dn 

lim — — — 

m^o d log M 



lim 

M->0 



2M 



X [1 + (V9Mo)" 



2;j1 



exp 



where in the above we have defined 



y/BT[n + S\ 



B 



1"o 



AnA s kl 
(2tt) 3 



oo , (47) 



(48) 



Hence, we see that the number density of haloes per log- 
arithmic interval diverges, ]imjf_>o Mn(M) oo, but 
that the multiplicity function of haloes asymptotes to a 
constant value, lirriM-s-o M 2 n{M) — ► const. 

We show this behaviour more schematically in the left 
panel of Fig. [3] (dot-dashed curves), where we compare 
CDM and WDM mass functions. For large halo masses 
M > 10 13 /i~ 1 Af Q , the models are virtually indistinguish- 
able, even for the case of a WDM particle with mass 
rawDM = 0.25keV. However, for smaller halo masses, 
whilst there is a significant suppression, they are still 
highly abundant. 

Recent results from numerical simulations show that 
there is a strong suppression in the numbers of haloes 
with M < Mfs [35(. Also, there is some debate as to 
whether haloes below the cut-off can form at all [23|, [34| . 
We shall take the view that PS is correct for M > Mf s , 
but that below this mass one finds virtually no haloes. 
We model this transition in an ad-hoc manner using an 
error function as a smoothed step. Thus the mass func- 
tion in WDM becomes: 



dn 



WDM 



dlogAf 



1 + erf 



logm (M/M 6 ) 



Clog M 



dn 



WDM 



dlogM 



(49) 

where <7i og m controls the logarithmic width of the step, 
and Mf s controls the location of the step. In this work, for 
simplicity we shall take a\ og M — > 0, and so the function 
acts like a standard Heaviside function at Mf s . We show 
the results of this modification in Fig. [3] (solid lines). For 
this example we have chosen a\ og M = 0.5. 

The exact amount of suppression in the WDM mass 
functions, relative to CDM, can be better quantified 
through inspecting the right panel of Fig. [3] Here we 
show the fractional difference between the WDM and 
CDM models versus halo mass-scaled in units of Mf s . 
The plot shows that there is roughly a ~ 50% suppres- 
sion in the abundance of haloes with M < 100 — 300 Mf s 
in all of the WDM models considered. 



B. Halo bias 

The halo model also requires us to specify how the 
centers of dark matter haloes of different masses cluster 
with respect to each other. As described by Eq. (|18p . for 
the Gaussian model, and at first order in the dark mat- 
ter density, halo and matter density perturbations can 
be related thr oug h a scale- independent bias factor b(M). 
Following [46J, [52J , an application of the peak-background 
split approximation enables one to calculate b(M) from 
a given mass function. For the ST mass function the 
Gaussian bias has the form: 



&St(^) = 1 



2p 



qv 2 — 1 

+ S c (z) [1 + [qu 2 Y] 



(50) 



where the parameters {p, q} are as in Eq. (j40]l. 

For the case of WDM, we shall assume that for haloes 
with M > Mf s the bias is well described by the same 
peak-background split model. Hence in model space: 
6 W dm(M) -> &cdm(M). 

In the left panel of Fig.@] we show how the bias of 
dark matter haloes in both CDM and WDM models 
depends on the halo mass. The figure shows the re- 
sults for the set of WDM particle masses: mwDM = 
{0.25, 0.5, 0.75, 1.0, 1.25} keV. As for the mass function, 
for high mass haloes with M > 10 13 /i _1 Mq, the bias is 
virtually indistinguishable from that obtained for CDM. 
For the case of lower masses, we see that the bias has 
a clear dependence on the mass of the WDM particle. 
For lower mass particles the bias of low mass haloes is 
significantly boosted with respect to CDM. 

The right panel of Fig. 2] shows this last point in more 
detail. Here we plot the ratio of the bias in the WDM 
model with respect to that of the CDM model, and we 
scale the halo mass in terms of the free-streaming mass- 
scale Mf s . We now see that there is roughly a 10% (2%) 
boost in the bias for haloes with M < 100 Mf s and with 
a WDM particle mass mwDM = 0.25 (1.25)keV. 

The left panel of Fig. 2] also shows how the bias 
of the smooth mass distribution varies with respect 
to wiwdm- This was calculated using Eq. ([35j). and 
assuming the ST mass function and bias parame- 
ter from above. We find that the smooth compo- 
nent is anti-biased with respect to the total matter, 
with b s e [0.730, 0.718, 0.716, 0.715, 0.715] for the cases 
mwDM = {0.25, . . . , 1.25} keV, respectively. Note that as 
the particle mass increases the bias of the smooth com- 
ponent saturates. This arises due to the fact that, for 
M/M* < 1, v -> and hence from Eq. (1501 . 



lim 6st(^) 
v— >o 



1 



[2p- 
D(z) 



1] /Sc(z) 
0.4 



1.686 



0.7 (z = 0) . (51) 
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FIG. 5: Density profiles of CDM and WDM haloes as a function of radius scaled in units of the characteristic scale radius r s . Top 
left, top right, bottom left and bottom right panels show the results for haloes with masses: M £ {10 15 , 10 13 , 10 11 , 10 9 }/i _1 Mq. 
In all panels: the black line denotes the CDM profile; the (red, green, blue, cyan, magenta) lines denote the WDM profiles for 
masses hiwdm G {0.25, 0.5, 0.75, 1.0, 1.25} keV, respectively. The solid lines denote the case where there are no thermal relic 
velocities; the dot-dash lines denote the case where profiles have been smoothed with a Gaussian filter of the scale of the mean 
relic velocity speed l V[ ~ v/Ho ~ 3.15vo/Hq. In each panel we also indicate the scale radius that the halo would have for the 
various WDM particle masses considered, where from top to bottom we show this for decreasing ttiwdm- 



C. Halo density profiles 

The density profiles of CDM haloes have been studied 
in great detail over the last two decades. A reasonably 
good approximation for the spherically averaged density 
profile in simulations is the NFW model [Hj]: 



P^Z(r\M)=pA c (M) y(l + y) 



y = r/r s , (52) 



where the two parameters are the scale radius r s and the 
characteristic density A c . If the halo mass is M v j r = 
47rr^ ir 200p/3, then owing to mass conservation there is 
only one free parameter, the concentration parameter 
c(M) = r vir /r s . Hence, 



A C (M) = (200/3)c 3 [log (1 



c/(l + c)]- 



(53) 



The concentration parameter can be obtained from the 
original model of NFW, but we shall instead use the 



model of Bullock et al. [53[. Note that for our halo model 
calculations we shall correct c(M) for the fact that the 
definitions of the virial mass in Bullock et al. and Sheth 
& Tormen are different [54]. 

In the case of the WDM model, there are two major 
effects which modify the density profiles when compared 
to CDM. Firstly, if one adopts either the NFW [|2] or 
Bullock et al. [53] approach for calculating the concen- 
tration parameters, then one finds, owing to the suppres- 
sion of small-scale power, that the collapse redshifts of 
the haloes are significantly affected. Thus a given halo 
of mass M in the WDM model, will be expected to col- 
lapse at a later time than in the CDM model. Since halo 
core density is related to the density of the Universe at 
the collapse time, then one expects that core densities in 
WDM models will be suppressed, and haloes on average 
to be less concentrated. 

We explicitly demonstrate this point in Fig. [5] 
The top left, top right, bottom left, and bottom 
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right panels show density profiles for haloes of mass 
Me [10 15 10 13 , 10 11 , 10 9 ] /i -1 M , respectively. The 
CDM case is shown as the black solid line, and the case 
where the mass of the dark matter particle is decreased, 
is shown as the series of curves with decreasing ampli- 
tude. Clearly, smaller mass haloes are more significantly 
affected by the absence of small-scale power. However, 
we do note that higher mass haloes are not immune to 
this, there being signs of density suppression in haloes 
with M ~ 10 13 h~ l M Q for m WDM ~ IkeV. 

Secondly, it has been argued, owing to the preservation 
of Liouville's theorem and the existence of a relic thermal 
velocity distribution for the dark matter particles, that 
there is a fundamental limit on the size of core densities 
in phase space [H, EH- If this is correct, then cuspy den- 
sity profiles are prohibited in the WDM model. In order 
to mimic the effect of the phase space constraint on the 
profiles, we take a qualitative approach and simply con- 
volve the NFW profile (obtained from computing collapse 
times with the WDM power spectrum) with a Gaussian 
filter function of radius the mean thermal velocity scale 
of the WDM particles. 

This can be calculated as follows: Let us suppose that 
the WDM particles are leptons and so obey the Fermi- 
Dirac (FD) statistics. Also, we assume that the gas de- 
coupled from the rest of matter whilst still relativistic, 
i.e. the energy of a given particle is E ~ cp. Thus its 
phase space distribution remains unchanged, except cool- 
ing adiabatically through the expansion of the Universe. 
Hence, the FD distribution function can be written: 



F(v) oc (1 + exp[t>/uo]) 



(54) 



Fitting to the results from the evolution of Einstein- 
Boltzmann codes gives the characteristic velocity (tem- 
perature) to be [23l | : 



v (z) w 0.012(1 + z) 



^WDM 


1/3 


h 


0.3 




0.65 



2/3 



"ffWDM" 


-1/3 


' m WDM' 


I 1.5 - 




. IkeV J 



-4/3 



kms 1 ,(55) 



where c/wdm is the number of degrees of freedom per 
particle, and we take this to be similar to that of a 
light neutrino species, hence jwdm ~ L5 (23[. Next, 
in order to calculate statistics, we also require the den- 
sity of momentum states, for a quantum gas this is: 
N(p) oc d 3 p/h 3 oc d 3 v oc N(v). Hence, the mean and 
rms velocities can now be calculated: 



TABLE I: Variation of the present day (z = 0) relic velocity 
statistics with WDM particle mass. Columns are: particle 
mass, characteristic velocity of FD gas, mean velocity, rms 
velocity and the smoothing scale for profiles. 



flWDM 

[keV] 


vo 
[kms -1 ] 


«7 
[kms -1 ] 


[kms -1 ] 


[/i _1 kpc] 


0.25 


0.065 


0.204 


0.234 


2.04 


0.50 


0.026 


0.082 


0.094 


0.82 


0.75 


0.015 


0.047 


0.054 


0.47 


1.00 


0.010 


0.032 


0.036 


0.32 


1.25 


0.008 


0.025 


0.029 


0.25 



C(3) fa 1.202 and ((5) « 1.037 are Riemann Zeta func- 
tions. Evaluation of Eqs ([56]) and (|57|) gives: = 3.15t>o 

and yu^ = 3.60wo- Table Q] presents the velocity statis- 
tics for our fiducial set of WDM particle masses. 

Finally, the density profile smoothed on the scale of 
the mean relic velocity field, is given by 



«(r|M) 



d 3 k 



where L 



(2tt) 3 

W T (z)/H Q , W G (k\l Vr ) 



MU^(k\M)W G (k\l Vi )j a (kr) 



exp 



(59) 

-{kl v y/2\ and 



where L^dm (r\M) = p c ^$ [ (r\M)/M is the mass normal- 
ized profile. We truncate the CDM profile at the virial 
radius and for the NFW model the Fourier transform can 
be written as [43| : 

tf~\TTNFVfru\ n*\ sin (fccr s ) 

f(c)U CDM (k\M) = - M1 + c) 

+ cos [kr s (l + c)] {Q [kr s (l + c)\ - C 4 [fcr s ]} 
+ sin [fcr„(l + c)] {S % [Ar.(l + c)] - &[fcr 8 ]}(60) 



where /(c) = [log (1 + c) — c/(l + c)] and where Si and 
Cj are the standard sine and cosine integrals. 

In Fig. [5] we show the impact of relic velocities on the 
density profiles. If the relic velocities are present and act 
to simply smooth the density field, then, by fiat, we see 
that halo cusps are suppressed. Importantly, the figure 
shows, for the range of WDM particle masses considered, 
that only haloes with M < 10 11 show any noticeable 
effects due to the relic velocities on scales r > O.lr,. 



RESULTS 



= A 



A 



dvvN(v)F(v) = 



dvv 2 N{v)F{v) 



7tt 4 



180C(3) 
1505) 
C(3) 



(56) 
(57) 



In order to obtain the last equalities in both Eqs ([56 
and (1571). we have used the normalization condition, 



1 

A 



dvv' 



1 = §C(3)w, 



(58) 



A. Mass power spectra 

In Fig. [6] we compare the WDM and CDM nonlinear 
matter power spectra obtained from the halo model cal- 
culations. The figure shows the contribution of each 
of the halo model terms to the WDM spectrum. For 
the purposes of clarity we show this for one case only, 
a WDM particle with m^DM = 0.25 keV. This is in- 
structive, as it shows that on large scales the power is 
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FIG. 6: Comparison of the nonlinear matter power spec- 
trum in the halo model for a particular WDM model 
(towdm = rn x — 0.25keV) and the CDM model, as a func- 
tion of wavenumber. The solid black and red lines denote the 
total CDM and total WDM power spectra in the halo model. 
The dashed, dot-dash, dotted and triple dot-dashed lines de- 
note, P BS , P B h, Pm, and P2H, respectively. The solid blue and 
red lines correspond to results computed with and without 
the relic velocities of the WDM particles - note that these 
lines are indistinguishable in the plot. 



totally determined by the weighted contributions P ss , 
P s h and -P2H- However, on small scales the power is en- 
tirely dominated by the haloes that exist above the free- 
streaming mass-scale. The integrand of the 1-Halo term 
is n(M)M 2 \U(k\M)\ 2 . Since there are virtually no differ- 
ences between the mass function in WDM and CDM for 
haloes above M ~ 1O 13 /i _1 M0, the apparent suppres- 
sion in power comes entirely from the suppression in the 
abundance of small haloes and the change in the profiles. 

In Fig. |6] we also show the importance of the relic ve- 
locities on the matter power spectrum. As can clearly be 
seen, at present times the halo model power spectra with 
and without the relic velocity effects are indistinguish- 
able for scales k < lOO/iMpc" 1 . On these scales, the 
modifications to the power spectrum due to uncertainties 
in the baryonic physics are significantly more important. 
We therefore assume that the present day effects of relic 
velocities will be unimportant for constraining the WDM 
particle from large-scale structure probes. 

In the left panel of Fig. we show the dependence of 
the nonlinear matter power spectra on the WDM particle 
mass. We see that, whilst the linear theory power spectra 
are significantly damped, the nonlinear halo model pre- 
dictions show that a large amount of small scale power 
is generated. As noted above this arises due to the col- 
lapse of haloes above the free-streaming mass-scale. The 



bottom panel of Fig. [7] allows us to quantify the differ- 
ences in greater detail. Here we show the ratio of the 
linear WDM and CDM spectra, and also the ratio of the 
nonlinear WDM and CDM spectra. From this we under- 
stand that the relative difference in the nonlinear model 
is significantly less than the difference between the linear 
theory predictions. For our cosmological model, the ef- 
fect of going from the CDM model to WDM model with 
particle mass towdm ^ IkeV, is to cause a ~ 10% sup- 
pression of power at k > 30 h Mpc -1 . 

In the right panel of Fig. we show the nonlinear pre- 
dictions for the WDM power spectra from the halof it 
code of [56j . This semi-empirical method was designed to 
match CDM power spectra. It is not clear whether this 
method can be accurately extrapolated to model nonlin- 
ear WDM spectra. Comparing the left and right pan- 
els of the figure, we see that the halo model predicts a 
stronger suppression of nonlinear power than halof it, 
and also a greater difference between the predictions for 
each WDM particle. It will require hi-resolution iV-body 
simulations, to disentangle which of these two models 
provides the better description for the power spectrum 
in WDM models. 

At this juncture, one might question the validity of 
the halo model predictions given our assumptions con- 
cerning the cut-off mass scale M cut = Mf s and also our 
adoption of the NFW profile for these same haloes. In 
Appendix IA II and IA 21 we show that the predictions for 
the power spectrum are insensitive to the precise value of 
Mcut (no changes are found over two-orders of magnitude 
in mass around Mf s ), and the exact shape of the profiles 
of these objects. More quantitatively, the predictions are 
modified by < 1% for scales k < lOO/iMpc^ 1 and for 
WDM particle masses mwDM > 0.25 keV. 



B. Weak lensing convergence power spectra 

We now turn to the question: Can the modifications 
to the dark matter model, of the kind considered in this 
paper, be detected with next generation surveys? To 
answer this, we shall focus on weak gravitational lensing 
by cosmic LSS as the most direct observable. 

Weak lensing is the magnification and shearing of 
light rays from distant source galaxies as they propagate 
through the inhomogeneous matter in the Universe on 
their way towards the observer. Owing to our ignorance 
of the intrinsic shapes of the unlensed galaxies, these dis- 
tortions can only be measured through the correlations 
in the shapes of the source galaxies with each other. Fol- 
lowing [57] . the weak lensing angular power spectrum of 
convergence is given by 



C, 



4c 1 



XH W (y) 
a (X) 



I 



Da(x)' 



where x is the comoving geodesic distance from the ob- 
server to redshift z, and Da(x) is the comoving angular 
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FIG. 7: Comparison of the halo model predictions for the WDM and CDM matter power spectra as a function of wavenumber. 
Left and right panels: results as obtained from the Halo Model and halof it, respectively. Top panels: absolute power. Solid and 
dashed lines denote the nonlinear and linear theory predictions, respectively. Black lines denote CDM; and colors {red, green, 
blue, cyan, magenta} denote WDM particle masses rawDM = mx £ {0.25, 0.5, 0.75, 1.0, 1.25} keV. Bottom panels: ratio of 
WDM and CDM power spectra. Line styles unchanged. 



diameter distance to redshift z. The survey weight func- 
tion is defined: 



W( X ) = / dx'G(x') 



Da(x'-x) 
Da(x') 



(62) 



The function G(x) is the normalized source distance 
distribution function, which tells us the number of 
sources per unit distance interval dx- For sim- 
plicity we shall assume the source distribution is a 
Dirac delta function at a single source plane, hence 
G(x) — 3° ix ~ Xs)- In this case the survey weights sim- 
plify to W(x) = D A (x s - x)/Da{Xs) and we have, 



Ci 



4c 4 



dx 



Da(Xs-x) 
Da(x s ) 



P 



I 



Da(xY 



Direct numerical integration of the above expression can 
be very slow, since it requires the evaluation of the halo 
model integrals for many points in the (k, a) plane, for a 
given I. The most efficient way to overcome this problem, 
is to generate a bi-cubic spline [HH fit to P(k, a) evaluated 
on a cubical grid of log k and log a. We found a mesh of 
size 100 x 50 gave good results, with k <E [0.005, 100.0] 
and a € [0.5, 1.0]. Once the bicubic spline is generated, 



the evaluation of the halo model C/ takes of order the 
same time to compute as the linear C; calculation. 

We now make a rough assessment of the expected er- 
rors on the convergence power spectrum for a future weak 
lensing survey. To do this, we assume that the spherical 
harmonic multipole coefficients are independent Gaus- 
sianly distributed variables. In this case the covariance 
matrix for a measurement of the power spectrum in band 
powers I and I' can be written 471 ] : 



Cav[Ci,C v ] = 



1 



1 



/sky N(l m , l n ) 



-1 2 



°Ll' 



(64) 



where / s ky is the fraction of sky covered, n is the angular 
number density of galaxies on the sky, cr 7 is the shape 
noise error, which is the error on a given measurement 
of the ellipticity of a galaxy. N(l mi l n ) is the number 
of independent multipoles used to generate the estimate 
of C'i in a given band power. For a given band power 
I G [l m , l n ], we find that 



N(l m ) In) 



E 

l = lrr, 



21 + 1 _ (l n + l m + 1) (l n -l m + l) 



(65) 
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FIG. 8: Comparison of the weak lensing convergence power spectrum in the CDM and WDM models, as a function of angular 
multipole for source galaxies at z s = 1. Left and right panels: predictions obtained using the halo model and halofit, 
respectively. Top panels: Absolute power. Dashed and solid lines show the predictions from linear theory and from the 
nonlinear halo model. Black lines represent results for CDM. Colored lines with {red, green, blue, cyan, magenta} denote 
results for the WDM model with particle masses ttiwdm = mwDM £ {0.25, 0.5, 0.75, 1.0, 1.25} keV. Bottom panel: Ratio of 
the WDM model predictions to the CDM predictions. Lines styles and colors are as above. 



For our fiducial survey we assume future full-sky weak 
lensing mission, like the proposed Euclid [59j or 
LSST [60| surveys: / s k y = 0.5, n — 35 galaxies per square 
arc minute, <r 1 — 0.35 and we take 50 bins in log-Z in the 
range (1 < / < 10 5 ). 

The left panel of Fig. [5] shows the convergence matter 
power spectrum for the linear and nonlinear halo model 
predictions for both CDM and the WDM models. For 
large angular scales, / < 10 2 , the linear and nonlinear 
predictions agree and there is no difference between CDM 
and WDM. However, on smaller scales we see significant 
departures between the predicted spectra. The impor- 
tance of the nonlinear corrections is very apparent: by 
I ~ 10 3 , there is an order of magnitude difference between 
the linear and nonlinear predictions, and by I ~ 10 4 , the 
difference is more than two orders of magnitude for the 
WDM models. As expected, we also see that the mea- 
surement errors are largest on very large and small scales, 
due to cosmic variance and shape noise, respectively. 

The bottom panel of Fig. [8] shows the relative differ- 
ences between the WDM and CDM models. We see 
clearly that nonlinear evolution reduces the differences 
between the WDM models and CDM. However, the error 



bars show that the WDM models and the CDM should 
be distinguishable at high significance. 

The right-hand panel of Fig. [5] shows the same but for 
the case where the nonlinear power is computed using 
halofit. The main differences are: for / ~ 10 3 the pre- 
dictions from both methods are roughly comparable, but 
for I > 10 4 the halo model has at least a factor 2 less 
signal than that obtained from halofit. The implica- 
tions of this are that if the halofit is correct and the 
halo model is wrong, then this will diminish the ability of 
future weak lensing surveys to constrain the mass of the 
WDM particle. Nevertheless, as was shown recently in 
(39l ] , reasonably good constraints may still be achievable. 

Note that for these 'signal-to-noise' predictions, we 
have assumed that all of the source galaxies lie at z s = 1. 
In the next section we perform a more complete Fisher 
matrix analysis with a full redshift distribution. 



C. Fisher matrix forecast for ttiwdm 

Finally, we make a forecast for the lower limit on 
mwDM obtainable with our fiducial (Euclid- like/LSST- 
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FIG. 9: We plot 65% and 95% likelihood contours, marginalised over all other parameters, for the fiducial (Euclid/LSST-like) 
survey (green) as well as the fiducial survey combined with Planck (blue). Our underlying fiducial model is the best fit WMAP7 
ACDM model. 



like) survey. Owing to these future missions including 
extensive photometric and spectroscopic redshift survey 
components, we now take advantage of this expected 
knowledge and assess the performance of a tomographic 
lensing analysis [see for example [6l[ . 

We adapt the analysis of the previous section as fol- 
lows. We assume an extended source redshift distribution 
for G(z) in Eq. , and for this we take the same form 
as in [30]. We then divide the distribution into 10 sepa- 
rate tomographic redshift bins. The lensing convergence 
power spectrum from the cross-correlation of sources in 
bins i and j is then given by: 



a 2 (x) 



Da(x) ,a J ' 

(66) 



where the weight function Wi(\) has the form: 



Wi( X ) 



= / ftoM^/mai(x,Xi) ^XsPs(z) £^ X ^ X {\ < Xi+l) 

\0 (x>Xi+i) 

(67) 

and where Wq = 3/2f2 TTli o^o an d Ps{z) is the source 
redshift distribution. Assuming that the underlying den- 
sity field obeys Gaussian statistics, then it can be shown 
that measurements of the convergence power spectra for 
different tomographic bins are correlated provided they 



share the same spherical harmonic multipole I. In this 
case the covariance matrix can be written [471: 



Cov 



1 



1 



/sky ^ (^mi ^n) 



°Ll' 



obs / 



^fobs . 



where we defined C°^(l) = C^j){l) + S^o^/m. In the 
above expression we have again ignored contributions to 
the covariance from higher order spectra [62]. 

Assuming that the likelihood function for obtaining 
measurements of the convergence power spectrum in sev- 
eral band powers and for several tomographic bins is well 
described by a multivariate Gaussian distribution with 
covariance matrix as given by Eq. (|68[) . then the Fisher 
matrix can be written [471: 



T ^^l^l^ Cov 



I X,Y 



C ( x)(l) 
dp 

(69) 

where a and (3 are elements of the vector of cosmologi- 
cal model parameters upon which our measurements de- 
pend. X and Y label the distinct pairs of tomographic 
bins, i.e. 16 {(1,1),..., (1,10), (2, 2),..., (2, 10),...}. 
The marginalized variance and covariances of the cosmo- 
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logical parameters can be obtained in the usual way [for 
further details see |(0|: o 2 aa = and o 2 afi = T~^. 

For our lensing analysis we choose to vary five parame- 
ters p € Wvvdm' ^m.Oi n s, A s , T}: the inverse WDM par- 
ticle mass; the total matter density today; the primordial 
perturbation spectral index; the primordial spectral am- 
plitude; and the matter power spectrum shape parame- 
ter, r = e _anb fl m fo. Note that, since our fiducial model 
is ACDM with fiducial parameters given by the WMAP7 
data [sij , we choose to work with Wwdm as a parameter 
in order to avoid setting the fiducial model parameter 
to infinity. However, redefining this variable does not 
avoid the strong variation of the error on m-yvDM with 
the chosen fiducial model. Furthermore, and most prob- 
lematically for such an analysis, the likelihood function is 
fiat with respect to the variation of this parameter around 
the maximum likelihood point. For this reason we fit half 
a Gaussian curve on top of this likelihood function and 
find single tailed errors as described in Markovic et al. 

Figure |H] presents the 1- and 2-D likelihood contours 
centered on our fiducial model. For the 2-D likelihood 
contours the 1- and 2-a contours are denoted by the light 
and dark green shaded regions, respectively. 

We then combine our lensing forecast with Planck 
65] priors, i.e J"^ tal = J"^ nsing + •F^ anck . We use the 
Planck Fisher matrix for the fiducial parameter set pre- 
scribed by [66]. We marginalize over the optical depth 
and the baryon density and perform a parameter trans- 
formation as in (oTI . l68j . As in [39| , we assume that we 
will only have the 143 GHz channel for CMB data. Fur- 
thermore, we assume that WDM with particle masses of 
the order of a keV, does not leave an observable signature 
in the CMB and hence set the corresponding entries in 
the Planck Fisher matrix to zero. The result of adding 
the Planck priors can be seen in the Fig. [H] as the blue 
contours and lines. 

Finally our forecast for the minimum limit on the 
WDM particle mass, assuming that the CDM scenario 
is correct, from a Euclid/LSST-like survey is wwdm > 
1.4keV. On combining our fiducial lensing survey and 
Planck data, the constraint tightens to m-yvDM > 2.6keV. 
We note that this limit is comparable to that found in 
39]. However, our analysis is more conservative since 
we consider only multipoles I < 5000, compared to 
I < 20,000 used in the study by @. Moreover, at the 
median redshift of the Euclid survey, z m = 0.9, I < 5000 
corresponds to approximately fc max = 1.6ft.Mpc _1 . The 
main reasons for excluding small scales from the analy- 
sis are: one, to avoid the range of scales where unknown 
baryonic physics effects inside clusters may be a signifi- 
cant source of bias for interpreting the data [Hj]; two, the 
assumption of a decorrelated covariance matrix for differ- 
ent Us breaks down due to nonlinear mode coupling in the 
density field, resulting in a loss of information ffuL FrO ] . 



VI. CONCLUSIONS 

In this paper we have explored a new approach to mod- 
elling the nonlinear evolution of two-point mass cluster- 
ing statistics in the WDM scenario. The model is based 
on the phenomenological halo model approach. 

In SjlT]we gave a brief overview of how linear structure 
formation is changed in the WDM scenario. We defined 
a free-streaming mass-scale Mf s , below which halo for- 
mation is likely to be strongly suppressed. 

In £11111 we developed the halo model of large-scale 
structure. We extended the model by including two com- 
ponents of clustered mass in the Universe. That which is 
in collapsed dark matter haloes and that which is in the 
form of smooth uncollapsed mass. This extended model 
requires that we understand three correlation functions: 
the halo-halo, smooth-halo and smooth-smooth correla- 
tion functions. The halo-halo clustering can be calculated 
as in the standard halo model. The clustering of the 
smooth matter is biased with respect to the total mass. 
We showed that, provided one understands the halo-halo 
correlation function, then the bias of the smooth compo- 
nent can be derived from simple conservation arguments. 

In ijlVl we explored how the standard halo model ingre- 
dients, i.e. halo mass function, halo bias and density pro- 
files are modified in the context of the WDM model. We 
showed that an application of the standard mass function 
theory, using the WDM power spectrum, leads to a 50% 
suppression in halo abundance relative to CDM for mass- 
scales M ~ lOOMfs. For our halo model calculations we 
assumed that no haloes formed below M{ s . 

For the halo bias, we assumed that the standard 
peak-background split calculations applied equally to the 
WDM model. Under this assumption, we showed that 
the bias of haloes with M < 100Mf s can be larger by 
~ 20% for WDM models than CDM ones, for m WDM = 
0.25 keV. Otherwise, halo bias is not affected. We also 
computed the bias for the smooth component and found 
that it was anti-biased with respect to the total mass, 
typically with a bias of order b ~ 0.7, for the range of 
WDM particle masses that we considered. 

In modelling the impact of WDM on the density pro- 
files, we took into account two effects: firstly, the change 
in core density and concentration due to the lower forma- 
tion redshift in WDM models; secondly, the suppression 
of density cusps due to relic thermal velocities of the 
WDM particles. We found that the former effect is the 
most important. Relic velocities only affect the density 
structure on scales r < 1 ft kpc, these are irrelevant for 
weak lensing analysis and so can be neglected. 

In $V]we computed the nonlinear matter power spec- 
trum in the halo model and showed that, whilst the linear 
power is strongly suppressed on small scales, nonlinear 
evolution generates significant small-scale power. Rela- 
tive to CDM the nonlinear WDM spectra are suppressed 
by 10% on scales k ~ 1 h Mpc -1 for a WDM particle with 
mass mwDM = 0.25 keV. For a heavier WDM particle 
the suppression is much less, for mwDM = 1-0 keV the 
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suppression is only ~ 10% on scales k ~ 20/iMpc -1 . We 
compared the predictions from the halo model with those 
from the halof it prescription (56|. We found that the 
latter showed more small scale power, but that the differ- 
ences between results for different WDM particle masses 
was not as apparent as in the case for the Halo Model. It 
will require high resolution iV-body simulations to vali- 
date which of the two models is more realistic. 

We then computed the impact of a set of WDM parti- 
cle masses on the expected convergence power spectrum 
from a future weak lensing survey. We also calculated the 
expected errors, assuming that each spherical harmonic 
is an independent Gaussian variable. We found that, if 
our halo model is correct, then it would be possible to 
differentiate between different WDM candidates with a 
systematics free, full-sky, weak lensing survey like Euclid 
III or LSST [eg. 

Finally, using the Fisher matrix approach we then fore- 
casted the lower limit on the WDM particle mass that 
would obtain from a Euclid like weak-lensing mission. 
We made the conservative assumption that multipoles 
I < 5000 can only be safely interpreted, and we found 
the lower limit towdm > 1.4 keV for Euclid only. On 
combining this with a forecast of the Planck CMB data, 
the constraint tightened to towdm > 2.6 keV. Our find- 



ings are in agreement with the earlier study of [39] , how- 
ever these earlier results were obtained assuming the less 
conservative multipole range I < 20, 000. 

CDM and WDM are equally plausible candidates for 
dark matter. We should therefore invest resources in to 
future cosmological experiments that allow us to place 
constraints on the mass of the dark particle. We expect 
that a full-sky weak lensing survey combined with the 
Planck data will greatly help achieve this objective. 
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Appendix A: Testing the assumptions 

The halo model predictions hinge upon a number of 
assumptions. We shall now test some of these. 



1. Impact of Cut-off mass scale on predictions 

Perhaps, the most questionable assumption is that we 
assume only haloes with M > Mf s exist and contribute 
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FIG. 10: Dependence of the halo model predictions for 
the WDM power spectra on the minimum mass scale 
A/cut, as a function of wavenumber. Top panel: ab- 
solute power. Solid lines of increasing thickness de- 
note the predictions with M cut G {0.1, 1.0, 10.0}A/ fs , re- 
spectively. Black lines denote CDM; and colors {red, 
green, blue, cyan, magenta} denote WDM particle masses 
mwDM = rnx 6 {0.25, 0.5, 0.75, 1.0, 1.25} keV. Bottom pan- 
els: ratio of WDM power spectra with varying cut-off mass 
scales to the fiducial WDM predictions with M cut = Mf s . 



to the clustering strength on small scales. As was dis- 
cussed earlier in $n] and PV A| there are good reasons 
for supposing that this is qualitatively correct. However, 
the exact relation must be determined from numerical 
simulations. Nevertheless, we may explore the impor- 
tance of this assumption for the predictions using our 
phenomenological model. 

In the top panel of Figure HD1 we show how the nonlin- 
ear power spectra for our standard set of WDM models 
varies with the cut-off mass scale. We consider three 
values for M cut € {0.1, 1.0, 10.0}Mf S) and in the figure 



these are depicted as the continuous lines of the same 
colour but with increasing thickness. As can be seen, no 
differences are apparent for the log-log plot. 

The bottom panel of Fig. [10] shows the results in 
greater detail, and here we present the fractional differ- 
ences between the power spectra with varying M cut and 
our fiducial case where M cut = Mf s . We find, for all of the 
scales and WDM models considered, that the differences 
are < 1%. The one exception is the ttiwdm = 0.25 keV 
model, for which we find that the differences deviate by 
> 1% on scales k > 30/iMpc _1 . 

Owing to the fact that these tests have covered two 
orders of magnitude in mass around the free-streaming 
mass scale, we are therefore satisfied that the halo model 
predictions are robust to changes in the fiducial value of 
M cut adopted in this work. 



Importance of density profile shape for haloes at 
the cut-off mass scale 



A further assumption that we have made is that all 
WDM haloes possess NFW density profiles and that the 
characteristic density and hence concentration parame- 
ters for these objects can be calculated as is done for 
CDM. However, it might be argued that haloes with 
M > M cut , have profiles which are structurally sub- 
stantially different from the NFW model. Of course, hi- 
resolution simulations are needed to provide the defini- 
tive answer to this question. However, we may use our 
phenomenological model to explore the importance of 
this assumption on the predictions. 

In fact, the analysis of the previous sub-section also 
helps clarify matters here too. Let us suppose that the 
NFW model provides an accurate description of haloes 
with M > M {s but for haloes with M > M fs that 
the profile is modified (32j. Let us consider the extreme 
case where the change to the halo profile in WDM is 
so violent that the halo effectively has no profile. This 
might simply be affected by removing the halo altogether 
from the calculation. As we have shown in the previous 
sub-section, if we add/remove haloes of mass one order 
of magnitude lower/higher than Mf s , the modifications 
to the predictions for the power spectrum are < 1% for 
wwdm > 0.25 keV. 

We thus conclude that the halo model predictions for 
the WDM power spectrum, on scales k < 100/iMpc -1 , 
are insensitive (< 1% changes) to the exact value of the 
cut-off mass scale M cut and the shapes of the density 
profiles near to Mf s , for towdm > 0.25 keV. 



